/*
Salil Benegal, Matt Motta
Data replication from 2020 ANES on determinants of Dunning Kruger attitudes wrt climate, vaccines, and covid
10/2/22 */

use "C:\Users\salil\Documents\2020 ANES\anes_timeseries_2020_stata_20210719.dta", clear
set more off

/* Preliminary data cleaning:
Partisanship, 7 and 3 point scales
Ideology
Racial resentment
Anti-Asian prejudice 
Sexism
Trust in experts
Education
Race
Sex
Income */

// Party ID
gen partyid3=.
replace partyid3=1 if V201228==1
replace partyid3=2 if V201228==2
replace partyid3=1 if V201230==3
replace partyid3=2 if V201230==1
replace partyid3=3 if V201230==2
label define partyid3 1 "Democratic" 2 "Republican" 3 "Independent/No pref"
label values partyid3 partyid3

tab V201228 partyid3
tab V201230 partyid3

// Ideology: this is problematic. 

tab V201200
gen ideology=V201200
recode ideology (-9=.m) (-8=.m) (99=.d)
tab ideology, m

recode ideology (2=1) (3=1) (4=2) (5=3) (6=3) (7=3)
tab V201200 ideology, m

label define ideology 1 "Liberal" 2 "Moderate" 3 "Conservative" 
label values ideology ideology

// Race
tab V201549x
gen race=V201549x
recode race(-9=.r) (-8=.d)
label define race 1 "White" 2 "Black" 3 "Hispanic" 4 "Asian/AAPI" 5 "Native American" 6 "Multiple/non-Hispanic"
label values race race

tab ideology if race==1, m
tab ideology if race!=1, m

// Sex
gen male=.
replace male=1 if V201600==1
replace male=0 if V201600==2
label define male 0 "Female" 1 "Male" 
label values male male

// Education
tab V201510, m 
gen educlevel=.
label define educ 1 "High school or less" 2 "Some college" 3 "Associates degree" 4 "Bachelors degree" 5 "Grad or higher", replace
replace educ=1 if V201510==1
replace educ=1 if V201510==2
replace educ=2 if V201510==3
replace educ=3 if V201510==4
replace educ=3 if V201510==5
replace educ=4 if V201510==6
replace educ=5 if V201510==7
replace educ=5 if V201510==8
label values educlevel educ
tab educlevel V201510, m
tab educlevel, gen(educ)

// age
gen age=V201507x
tab age, m
recode age(-9=.r)

gen agedem=.
label define agedem 1 "18-30" 2 "31-40" 3 "41-50" 4 "51-64" 5 "65+"
replace agedem=1 if age>17 & age<31 & !missing(age)
replace agedem=2 if age>30 & age<41 & !missing(age)
replace agedem=3 if age>40 & age<51 & !missing(age)
replace agedem=4 if age>50 & age<65 & !missing(age) 
replace agedem=5 if age>64 & !missing(age)

label values agedem agedem
tab age agedem, m

// 4 item racial resentment
// V202300: favors, 1 most racist response - need to reverse
// V202301: generations - 5 most racist response
// V202302: deserve - 5 most racist response
// V202303: tryhard - 1 most racist - need to reverse

gen favors=.
replace favors=1 if V202300==5
replace favors=2 if V202300==4
replace favors=3 if V202300==3
replace favors=4 if V202300==2
replace favors=5 if V202300==1

gen generations=.
replace generations=1 if V202301==1
replace generations=2 if V202301==2
replace generations=3 if V202301==3
replace generations=4 if V202301==4
replace generations=5 if V202301==5

gen deserve=.
replace deserve=1 if V202302==1
replace deserve=2 if V202302==2
replace deserve=3 if V202302==3
replace deserve=4 if V202302==4
replace deserve=5 if V202302==5

gen tryhard=.
replace tryhard=1 if V202303==5
replace tryhard=2 if V202303==4
replace tryhard=3 if V202303==3
replace tryhard=4 if V202303==2
replace tryhard=5 if V202303==1

factor favors generations deserve tryhard, pcf

gen resentment=(favors+generations+deserve+tryhard)/4
gen resentfactor=(resentment-1)/4

// 2 item hostile sexism
// V201639: innocent remarks, 1 most sexist response - need to reverse
// V201640: gain control - 1 most sexist - need to reverse

gen innocent=.
replace innocent=1 if V201639==5
replace innocent=2 if V201639==4
replace innocent=3 if V201639==3
replace innocent=4 if V201639==2
replace innocent=5 if V201639==1

gen control=.
replace control=1 if V201640==5
replace control=2 if V201640==4
replace control=3 if V201640==3
replace control=4 if V201640==2
replace control=5 if V201640==1

factor innocent control

gen hostilesexism=(innocent+control)/2
gen hostilefactor=(hostilesexism-1)/4

// trust in experts
gen trustexperts=.
replace trustexperts=V202308x if V202308x>0
label define trustexperts 1 "Trust ordinary people much more" 2 "Trust ordinary somewhat more" 3 "Trust both the same" 4 "Trust experts somewhat more" 5 "Trust experts much more"
label values trustexperts trustexperts

tab trustexperts V202308x, m
label var trustexperts "Trust in experts"

gen trustgovt=.
label define trustgovt 1 "Never" 2 "Some of the time" 3 "About half the time" 4 "Most of the time" 5 "Always"
replace trustgovt=1 if V201233==5
replace trustgovt=2 if V201233==4
replace trustgovt=3 if V201233==3
replace trustgovt=4 if V201233==2
replace trustgovt=5 if V201233==1
label values trustgovt trustgovt

// Dependent variables: 

tab V202555
tab V202556

// climate dunning kruger - wrong answer *and* moderate confidence or higher

gen dkclimate=.
replace dkclimate=0 if V202555==1
replace dkclimate=0 if V202555==2 & V202556==1
replace dkclimate=0 if V202555==2 & V202556==2

replace dkclimate=1 if V202555==2 & V202556==3
replace dkclimate=1 if V202555==2 & V202556==4
replace dkclimate=1 if V202555==2 & V202556==5

tab dkclimate
// only 9% of the ANES sample both misinformed and overconfident on the climate change happening q

logit dkclimate b3.partyid3 hostilefactor resentfactor i.educlevel male


// vaccines cause autism dk

tab V202553
tab V202554
tab V202554 if V202553==1

gen dkautism=.
replace dkautism=0 if V202553==2
replace dkautism=0 if V202553==1 & V202554==1
replace dkautism=0 if V202553==1 & V202554==2

replace dkautism=1 if V202553==1 & V202554==3
replace dkautism=1 if V202553==1 & V202554==4
replace dkautism=1 if V202553==1 & V202554==5

tab dkaut

logit dkaut b3.partyid3 hostilefactor resentfactor i.educlevel male

// covid developed in a lab - 

tab V202557
tab V202558
tab V202558 if V202557==1

gen dkcovidlab=.
replace dkcovidlab=0 if V202557==2
replace dkcovidlab=0 if V202557==1 & V202558==1 
replace dkcovidlab=0 if V202557==1 & V202558==2

replace dkcovidlab=1 if V202557==1 & V202558==3 
replace dkcovidlab=1 if V202557==1 & V202558==4
replace dkcovidlab=1 if V202557==1 & V202558==5

tab dkcovidlab 
logit dkcovidlab b3.partyid3 hostilefactor resentfactor i.educlevel male

// hydroxychloroquine (HCQ) views

tab V202559
tab V202560 if V202559==1

gen dkhcq=.
replace dkhcq=0 if V202559==2
replace dkhcq=0 if V202559==1 & V202560==1
replace dkhcq=0 if V202559==1 & V202560==2

replace dkhcq=1 if V202559==1 & V202560==3
replace dkhcq=1 if V202559==1 & V202560==4
replace dkhcq=1 if V202559==1 & V202560==5

tab dkhcq

logit dkhcq b3.partyid3 hostilefactor resentfactor i.educlevel male

gen dkscale=(dkclimate+dkautism+dkcovidlab+dkhcq)

// MAIN ANALYSIS BEGINS HERE

// figs 1 and 2

cd "C:\Users\salil\Dropbox\Work in progress\Motta Benegal\Appendix"

** Figure 1: RR, controls, full sample
set more off
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Vaccines and autism") xtitle("Effect size") saving(autism2, replace)
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Climate change") xtitle("Effect size") saving(climate2, replace)
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Covid lab origin") xtitle("Effect size") saving(lab2, replace)
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Hydroxychloroquine") xtitle("Effect size") saving(hcq2, replace)

graph combine autism2.gph climate2.gph lab2.gph hcq2.gph, col(2)


** Figure 2: predicted probabilities!
set scheme plottig
eststo clear
set more off
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store autism_dems
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==2 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store autism_reps
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store climate_dems
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==2 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store climate_reps
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store lab_dems
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==2 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store lab_reps
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store hcq_dems
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==2 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store hcq_reps

coefplot (autism_dems, label (Democrats) ) (autism_reps, label (Republicans) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Vaccines cause autism") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "autism_prob.gph", replace
	
coefplot (climate_dems, label (Democrats) ) (climate_reps, label (Republicans) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("World temps have not changed") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "clim_prob.gph", replace
	
coefplot (lab_dems, label (Democrats) ) (lab_reps, label (Republicans) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Covid-19 created in a lab") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%" .6 "60%" .7 "70%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "lab_prob.gph", replace
	
coefplot (hcq_dems, label (Democrats) ) (hcq_reps, label (Republicans) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Hydroxychloroquine effective treatment") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "hcq_prob.gph", replace
	
graph combine autism_prob.gph clim_prob.gph lab_prob.gph hcq_prob.gph, col(2)	


** Alt Figure 2: predicted probabilities but across ideo!
set scheme plottig
eststo clear
set more off
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel male age if ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store autism_libs
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel male age if ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store autism_cons
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel male age if ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store climate_libs
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel male age if ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store climate_cons
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel male age if ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store lab_libs
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel male age if ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store lab_cons
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel male age if ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store hcq_libs
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel male age if ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store hcq_cons

coefplot (autism_libs, label (Liberals) ) (autism_cons, label (Conservatives) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Vaccines cause autism") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "autism_prob.gph", replace
	
coefplot (climate_libs, label (Liberals) ) (climate_cons, label (Conservatives) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("World temps have not changed") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "clim_prob.gph", replace
	
coefplot (lab_libs, label (Liberals) ) (lab_cons, label (Conservatives) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Covid-19 created in a lab") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%" .6 "60%" .7 "70%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "lab_prob.gph", replace
	
coefplot (hcq_libs, label (Liberals) ) (hcq_cons, label (Conservatives) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Hydroxychloroquine effective treatment") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "hcq_prob.gph", replace
	
graph combine autism_prob.gph clim_prob.gph lab_prob.gph hcq_prob.gph, col(2)	






** Alternate Figure 2: Including RR, but only white respondents
set more off
logit dkautism resentfactor b3.partyid3 i.educlevel trustexperts male age if race==1 [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Vaccines and autism") xtitle("Effect size") saving(autism2, replace)
logit dkclimate resentfactor b3.partyid3 i.educlevel trustexperts male age if race==1  [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Climate change") xtitle("Effect size") saving(climate2, replace)
logit dkcovidlab resentfactor b3.partyid3 i.educlevel trustexperts male age if race==1  [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Covid lab origin") xtitle("Effect size") saving(lab2, replace)
logit dkhcq resentfactor b3.partyid3 i.educlevel male trustexperts age if race==1  [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Hydroxychloroquine") xtitle("Effect size") saving(hcq2, replace)

graph combine autism2.gph climate2.gph lab2.gph hcq2.gph, col(2)

** Figure 3: including group-centric ethnocentrism alongside RR, only white respondents
set more off
logit dkautism resentfactor ethnofactor* b3.partyid3 i.educlevel trustexperts male age if race==1 [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Vaccines and autism") xtitle("Effect size") saving(autism3, replace)
logit dkclimate resentfactor ethnofactor* b3.partyid3 i.educlevel trustexperts male age if race==1  [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Climate change") xtitle("Effect size") saving(climate3, replace)
logit dkcovidlab resentfactor ethnofactor* b3.partyid3 i.educlevel trustexperts male age if race==1  [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Covid lab origin") xtitle("Effect size") saving(lab3, replace)
logit dkhcq resentfactor ethnofactor* b3.partyid3 i.educlevel trustexperts male age if race==1  [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Hydroxychloroquine") xtitle("Effect size") saving(hcq3, replace)

graph combine autism3.gph climate3.gph lab3.gph hcq3.gph, col(2)


set more off
logit dkautism resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Vaccines and autism") xtitle("Effect size") saving(autism2, replace)
logit dkclimate resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Climate change") xtitle("Effect size") saving(climate2, replace)
logit dkcovidlab resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Covid lab origin") xtitle("Effect size") saving(lab2, replace)
logit dkhcq resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Hydroxychloroquine") xtitle("Effect size") saving(hcq2, replace)

graph combine autism2.gph climate2.gph lab2.gph hcq2.gph, col(2)

** Figure 3: looking at confidence only among those who endorse misinfo

tab V202553 V202554
recode V202554 (-1=.m) (-5=.m) (-6=.m) (-7=.m)
tab V202554, m
recode V202556 (-1=.m) (-5=.m) (-6=.m) (-7=.m)
tab V202555 V202556, m
recode V202558 (-1=.m) (-5=.m) (-6=.m) (-7=.m)
tab V202557 V202558, m
recode V202560 (-1=.m) (-5=.m) (-6=.m) (-7=.m)
tab V202560, m
tab V202559 V202560

// autism: 
set more off
ologit V202554 resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b] if V202553==1
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Vaccines and autism") xtitle("Effect size") saving(autism4, replace)
// climate
ologit V202556 resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b] if V202555==2
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Climate change") xtitle("Effect size") saving(climate4, replace)
// lab
ologit V202558 resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b] if V202557==1
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Covid lab origin") xtitle("Effect size") saving(lab4, replace)
// hcq
ologit V202560 resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b] if V202559==1
coefplot, drop(_cons) xline(0) xlabel(-2 -1 0 1 2) title("Hydroxychloroquine") xtitle("Effect size") saving(hcq4, replace)
graph combine autism4.gph climate4.gph lab4.gph hcq4.gph, col(2)

// or as a table

set more off
eststo clear
eststo autism: ologit V202554 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b] if V202553==1
eststo climate: ologit V202556 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b] if V202555==2
eststo covidlab: ologit V202558 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b] if V202557==1
eststo hcq: ologit V202560 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b] if V202559==1
 
esttab autism climate covidlab hcq using Table3.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Correlates of increased confidence only among respondents endorsing misinformation" )
	 label nonumbers
	//addnote("Dependent variable")//
	compress replace;
 #delimit cr

********************************************************************
********************************************************************
**************** APPENDIX STUFF/ROBUSTNESS CHECKS ******************
********************************************************************
********************************************************************
 
cd "C:\Users\salil\Dropbox\Work in progress\Motta Benegal\Appendix"

** Figure 1: PID, education, not RR
set more off
eststo clear
eststo autism: logit dkautism b3.partyid3 i.educlevel male age [pw=V200010b]
eststo climate: logit dkclimate b3.partyid3 i.educlevel male age [pw=V200010b]
eststo covidlab: logit dkcovidlab b3.partyid3 i.educlevel male age [pw=V200010b]
eststo hcq: logit dkhcq b3.partyid3 i.educlevel male age [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable1.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Correlates of Dunning-Kruger views" ) mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr

** Figure 2: Including RR
set more off
eststo clear
eststo autism: logit dkautism resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
eststo climate: logit dkclimate resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
eststo covidlab: logit dkcovidlab resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
eststo hcq: logit dkhcq resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable2.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Relationship between racial resentment and Dunning-Kruger views") mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr

** Alt Figure 2: RR, but only white respondents
set more off
eststo clear
eststo autism: logit dkautism resentfactor b3.partyid3 i.educlevel male age if race==1 [pw=V200010b]
eststo climate: logit dkclimate resentfactor b3.partyid3 i.educlevel male age if race==1 [pw=V200010b]
eststo covidlab: logit dkcovidlab resentfactor b3.partyid3 i.educlevel male age if race==1 [pw=V200010b]
eststo hcq: logit dkhcq resentfactor b3.partyid3 i.educlevel male age if race==1 [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable2a.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Relationship between racial resentment and Dunning-Kruger views (white respondent subsample") mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr


** Alt Figure 2: RR, but exclude Black respondents
set more off
eststo clear
eststo autism: logit dkautism resentfactor b3.partyid3 i.educlevel male age if race!=2 [pw=V200010b]
eststo climate: logit dkclimate resentfactor b3.partyid3 i.educlevel male age if race!=2 [pw=V200010b]
eststo covidlab: logit dkcovidlab resentfactor b3.partyid3 i.educlevel male age if race!=2 [pw=V200010b]
eststo hcq: logit dkhcq resentfactor b3.partyid3 i.educlevel male age if race!=2 [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable2b.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Relationship between racial resentment and Dunning-Kruger views (omit Black respondents") mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr


** Figure 3: including group-centric ethnocentrism alongside RR, only white respondents
set more off
eststo clear
eststo autism: logit dkautism resentfactor ethnofactor* b3.partyid3 i.educlevel male age if race==1 [pw=V200010b]
eststo climate: logit dkclimate resentfactor ethnofactor* b3.partyid3 i.educlevel male age if race==1  [pw=V200010b]
eststo covidlab: logit dkcovidlab resentfactor ethnofactor* b3.partyid3 i.educlevel male age if race==1  [pw=V200010b]
eststo hcq: logit dkhcq resentfactor ethnofactor* b3.partyid3 i.educlevel male age if race==1  [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable3.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Racial resentment, ethnocentrism, and Dunning-Kruger views (white respondent subsample") mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr

 ** Figure/table 4: looking at confidence, only among those endorsing misinfo
set more off
eststo clear
eststo autism: ologit V202554 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b] if V202553==1
eststo climate: ologit V202556 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b] if V202555==2
eststo covidlab: ologit V202558 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b] if V202557==1
eststo hcq: ologit V202560 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b] if V202559==1

esttab autism climate covidlab hcq using AppTable4.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Correlates of increased confidence only among respondents endorsing misinformation" ) mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr
 
** Robustness checks **
 
** Appendix table 1 recalculated with ideology
set more off
eststo clear
eststo autism: logit dkautism b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
eststo climate: logit dkclimate b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
eststo covidlab: logit dkcovidlab b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
eststo hcq: logit dkhcq b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable5.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Correlates of Dunning-Kruger views" ) mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr

** Appendix table 2 recalculated with ideology
set more off
eststo clear
eststo autism: logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
eststo climate: logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
eststo covidlab: logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]
eststo hcq: logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable6.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Relationship between racial resentment and Dunning-Kruger views") mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr

** Appendix table 3 recalculated with ideology
set more off
eststo clear
eststo autism: logit dkautism resentfactor ethnofactor* b3.partyid3 b2.ideo i.educlevel male age if race==1 [pw=V200010b]
eststo climate: logit dkclimate resentfactor ethnofactor* b3.partyid3 b2.ideo i.educlevel male age if race==1  [pw=V200010b]
eststo covidlab: logit dkcovidlab resentfactor ethnofactor* b3.partyid3 b2.ideo i.educlevel male age if race==1  [pw=V200010b]
eststo hcq: logit dkhcq resentfactor ethnofactor* b3.partyid3 b2.ideo i.educlevel male age if race==1  [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable7.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Racial resentment, ethnocentrism, and Dunning-Kruger views (white respondent subsample") mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr

** Appendix table 4 recalculated with ideology: looking at confidence, only among those endorsing misinfo
set more off
eststo clear
eststo autism: ologit V202554 resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b] if V202553==1
eststo climate: ologit V202556 resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b] if V202555==2
eststo covidlab: ologit V202558 resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b] if V202557==1
eststo hcq: ologit V202560 resentfactor b3.partyid3 b2.ideo i.educlevel male age [pw=V200010b] if V202559==1

esttab autism climate covidlab hcq using AppTable8.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Correlates of increased confidence only among respondents endorsing misinformation" ) mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr
 
// recalculate models 1-3 using ordinal DVs instead of binary 
 
label define dkscale 0 "Correct" 1 "Misinfo, low confidence" 2 "Misinfo, moderate+ confidence"

// vaccines cause autism dk

tab V202553
tab V202554
tab V202554 if V202553==1

gen dkautism2=.
replace dkautism2=0 if V202553==2
replace dkautism2=1 if V202553==1 & V202554==1
replace dkautism2=1 if V202553==1 & V202554==2
replace dkautism2=2 if V202553==1 & V202554==3
replace dkautism2=2 if V202553==1 & V202554==4
replace dkautism2=2 if V202553==1 & V202554==5

label values dkautism2 dkscale

// climate
gen dkclimate2=.
replace dkclimate2=0 if V202555==1 
replace dkclimate2=1 if V202555==2 & V202556==1
replace dkclimate2=1 if V202555==2 & V202556==2
replace dkclimate2=2 if V202555==2 & V202556==3
replace dkclimate2=2 if V202555==2 & V202556==4
replace dkclimate2=2 if V202555==2 & V202556==5

label values dkclimate2 dkscale

// covid lab
gen dkcovidlab2=.
replace dkcovidlab2=0 if V202557==2
replace dkcovidlab2=1 if V202557==1 & V202558==1 
replace dkcovidlab2=1 if V202557==1 & V202558==2
replace dkcovidlab2=2 if V202557==1 & V202558==3 
replace dkcovidlab2=2 if V202557==1 & V202558==4
replace dkcovidlab2=2 if V202557==1 & V202558==5

tab dkcovidlab2
tab dkcovidlab2 V202557
label values dkcovidlab2 dkscale

// hydroxychloroquine (HCQ) views

tab V202559
tab V202560 if V202559==1

gen dkhcq2=.
replace dkhcq2=0 if V202559==2 
replace dkhcq2=1 if V202559==1 & V202560==1
replace dkhcq2=1 if V202559==1 & V202560==2
replace dkhcq2=2 if V202559==1 & V202560==3
replace dkhcq2=2 if V202559==1 & V202560==4
replace dkhcq2=2 if V202559==1 & V202560==5
label values dkhcq2 dkscale

tab dkhcq2 
 
// now to redo the tables 
 
** Figure 1: PID, education, not RR
set more off
eststo clear
eststo autism: ologit dkautism2 b3.partyid3 i.educlevel male age [pw=V200010b]
eststo climate: ologit dkclimate2 b3.partyid3 i.educlevel male age [pw=V200010b]
eststo covidlab: ologit dkcovidlab2 b3.partyid3 i.educlevel male age [pw=V200010b]
eststo hcq: ologit dkhcq2 b3.partyid3 i.educlevel male age [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable9.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Correlates of Dunning-Kruger views with ordinal DV" ) mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr

** Figure 2: Including RR
set more off
eststo clear
eststo autism: ologit dkautism2 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
eststo climate: ologit dkclimate2 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
eststo covidlab: ologit dkcovidlab2 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]
eststo hcq: ologit dkhcq2 resentfactor b3.partyid3 i.educlevel male age [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable10.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Relationship between racial resentment and Dunning-Kruger views with ordinal DV") mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr

** Figure 3: including group-centric ethnocentrism alongside RR, only white respondents
set more off
eststo clear
eststo autism: ologit dkautism2 resentfactor ethnofactor* b3.partyid3 i.educlevel male age if race==1 [pw=V200010b]
eststo climate: ologit dkclimate2 resentfactor ethnofactor* b3.partyid3 i.educlevel male age if race==1  [pw=V200010b]
eststo covidlab: ologit dkcovidlab2 resentfactor ethnofactor* b3.partyid3 i.educlevel male age if race==1  [pw=V200010b]
eststo hcq: ologit dkhcq2 resentfactor ethnofactor* b3.partyid3 i.educlevel male age if race==1  [pw=V200010b]

 esttab autism climate covidlab hcq using AppTable11.rtf, nogap se starlevels(* .05 ** .01 *** .001) r2(%9.3f) ///
	stats(N, labels("N. of cases")) ///
	title("Racial resentment, ethnocentrism, and Dunning-Kruger views with ordinal DV") mtitle("Autism" "Climate" "Covid lab origin" "HCQ") ///
	label nonumbers compress replace 
 #delimit cr
 

 
 
 /////////////// fig 2
 
 ** Figure 2: predicted probabilities!
set scheme plottig
eststo clear
set more off
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==1 & ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store autism_dems
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==2 & ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store autism_reps
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==1 & ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store climate_dems
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==2 & ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store climate_reps
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==1 & ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store lab_dems
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==2 & ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store lab_reps
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==1 & ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store hcq_dems
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if partyid3==2 & ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store hcq_reps

coefplot (autism_dems, label (Democrats) ) (autism_reps, label (Republicans) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Vaccines cause autism") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "autism_prob.gph", replace
	
coefplot (climate_dems, label (Democrats) ) (climate_reps, label (Republicans) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("World temps have not changed") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "clim_prob.gph", replace
	
coefplot (lab_dems, label (Democrats) ) (lab_reps, label (Republicans) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Covid-19 created in a lab") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%" .6 "60%" .7 "70%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "lab_prob.gph", replace
	
coefplot (hcq_dems, label (Democrats) ) (hcq_reps, label (Republicans) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Hydroxychloroquine effective treatment") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "hcq_prob.gph", replace
	
graph combine autism_prob.gph clim_prob.gph lab_prob.gph hcq_prob.gph, col(2)	


** NEW Figure 2: predicted probabilities but across ideo
set scheme plottig
eststo clear
set more off
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store autism_libs
logit dkautism resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store autism_cons
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store climate_libs
logit dkclimate resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store climate_cons
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store lab_libs
logit dkcovidlab resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store lab_cons
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if ideo==1 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store hcq_libs
logit dkhcq resentfactor b3.partyid3 b2.ideo i.educlevel trustexperts male age if ideo==3 [pw=V200010b]
	margins, at(resentfactor=(0(.125)1))   atmeans post 
	estimates store hcq_cons

coefplot (autism_libs, label (Liberals) ) (autism_cons, label (Conservatives) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Vaccines cause autism") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "autism_prob.gph", replace
	
coefplot (climate_libs, label (Liberals) ) (climate_cons, label (Conservatives) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("World temps have not changed") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "clim_prob.gph", replace
	
coefplot (lab_libs, label (Liberals) ) (lab_cons, label (Conservatives) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Covid-19 created in a lab") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%" .6 "60%" .7 "70%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "lab_prob.gph", replace
	
coefplot (hcq_libs, label (Liberals) ) (hcq_cons, label (Conservatives) ), bylabel(Racial resentment) ///
	ytitle ("% Agreeing") ///
	xlab(0 "Low" .5 "Medium" 1 "High" , labsize(small)) xtitle ("Racial resentment") ///
	title ("Hydroxychloroquine effective treatment") ///
	ylab(0 "0%" .1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%", labsize(vsmall)) ///
	at recast(line) lwidth(*2) ciopts(recast(rline) lc(gray) lp(-) lw(vthin)) ///
	legend (col(2) pos(6)) 
	graph save "hcq_prob.gph", replace
	
graph combine autism_prob.gph clim_prob.gph lab_prob.gph hcq_prob.gph, col(2)	

